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ABSTRACT 

Context. Using the atmosphere as a detector volume, Imaging Air Cherenkov Telescopes (IACTs) depend highly on the 
properties and the condition of the air mass above the telescope. On the Canary Island of La Palma, where the Major 
Atmospheric Gamma-ray Imaging Cherenkov telescope (MAGIC) is situated, the Saharan Air Layer (SAL) can cause 
strong atmospheric absorption affecting the data quality and resulting in a reduced gamma flux. 

Aims. To correlate IACT data with other measurements, e.g. long-term monitoring or Multi-Wavelength (MWL) studies, 
an accurate flux determination is mandatory. Therefore, a method to correct the data for the effect of the SAL is required. 
Methods. Three different measurements of the atmospheric absorption are performed on La Palma. From the determined 
transmission, a correction factor is calculated and applied to the MAGIC data. 

Results. The different transmission measurements from optical and IACT data provide comparable results. MAGIC data 
of PG 1553+113, taken during a MWL campaign in July 2006, were analyzed using the presented method, providing a 
corrected flux measurement for the study of the spectral energy distribution of the source. 
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1. Introduction 

The Earth's atmosphere plays an important role in all 
ground-based observations because it is traversed by light 
from all objects observed. Scattering and absorption of light 
, by the atmosphere is variable and reduces the flux mea- 
sured at the telescope. On the Canary Island of La Palma, 
very strong atmospheric absorption can occur, in case the 
meteorologic phenomenon Saharan Air Layer (SAL) takes 
place. A warm air mass, consisting of mineral dust, travels 
from the Sahara westward; this becomes the so-called SAL, 
when it is undercut by a cool and moist air layer from the 
' sea. Located between 1.5 km and 5.5 km above sea level 
, (a.s.l.), the SAL, also known as Calima, can extend above 
the North Atlantic Ocean a s far as the United Stat es and 
the western Caribbean Sea (jDunion & Velden|[200l . 

Variable absorption can cause apparent flux variations, 
and the absolute flux level is in all cases reduced. To to cor- 
rect for this reduction, the atmospheric absorption is mea- 
sured. On La Palma, three measurements are available. The 
Carlsberg Meridian Telescope (CMT) provides nightly val- 
ues of extinction. From the optical data of the KVA and also 
from the MAGIC data itself, the atmospheric transmission 
can be determined. The results of these measurements are 
discussed and compared for the nights between July 15 th 
2006 and July 28 th 2006. 

A correction method was developed for the MAGIC tele- 
scope: From the absorption measurements, a factor is cal- 
culated and used in the analysis to correct for the effect of 
SAL. 
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2. Feasibility Study 



Before the development of the method, the feasibility of a 
correction based on the transmission measurement had to 
be checked. After explaining the principle of the Imaging 
Air Cherenkov Technique, the effect of the SAL on the 
Cherenkov data is explained, discussing in particular the 
vertical shower development. 



2.1. Imaging Air Cherenkov Technique 



Particles enter the atmosphere, interact with nuclei 
and produce electromagnetic or particle cascades. The 
Cherenkov light emitted by these secondary particles tra- 
verses the remainder of the atmosphere and is measured by 
a camera equipped with photon detectors. The morphol- 
ogy of the recorded shower images is used to distinguish 
the images of gamma-induced showers from those initiated 
by other particles, such as protons. The refore, the shower 
images are parameterized (Hillas 1985), and cuts in the 
corresp onding distributions a re used to suppress the back- 
ground (|Riegel fc B retz 2003). The energy spectrum is cal- 
culated by comparing the parameters of the measured to 
those of the simulated data, for which the energy is known. 
The most important parameter to reconstruct the energy of 
a primary photon is SIZE, i.e. the total amount of light in 
the shower image. Apart from minor dependences on other 
properties, the energy is mainly proportional to the SIZE. 
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2.2. Effect of Atmospheric Absorption 

The atmospheric absorption reduces the number of mea- 
sured photons and therefore the SIZE of the shower im- 
ages. Consequently, the energy of the primary particle is 
underestimated when comparing the SIZE of affected data 
to that of showers simulated for an ideal atmosphere. For 
an accurate energy reconstruction, the shower images have 
to be corrected for the loss of light due to absorption, or 
the atmospheric absorption has to be included in the sim- 
ulation. As different amounts of atmospheric absorption in 
the single nights would require different sets of simulated 
data, a correction method for the data is preferred. 

2.3. Dependence on the Height in the Atmosphere 

Based on optical data, the values for transmission provided 
by CMT and KVA correspond to light passing the whole 
atmosphere. The Cherenkov light of a shower originates in a 
height between 20 km a.s.l. and the ground. While the first 
interaction occurs at about 20 km a.s.l., the shower maxi- 
mum is at around 10 km a.s.l. Figure [1] shows the number of 
produced Cherenkov photons (upper plot) and the fraction 
of produced light (lower plot) versus height. From the inte- 
gral fraction of light (lower plot), it can be concluded that, 
for primary particles with an energy below 500 GeV, more 
than 90 % of the shower light is produced above 5.5 km a.s.l. 
Consequently, less than 10 % of the overall emitted light is 
absorbed by a slightly lower amount, because it does not 
traverse the entire SAL. To a first-order approximation, this 
small effect can be neglected, implying that a factor calcu- 
lated from the measured atmospheric transmission can be 
used to correct the affected data. 

3. Measurements of the Atmospheric Absorption 

Absorption measurements are available from three differ- 
ent instruments at the Observatorio del Roque de los 
Muchachos (ORM). 

3.1. CMT- Extinction Measurements 

The CMT ()Evans et all l2002f) is one of the optical tele- 
scopes in the ORM. Since 1984, it provides a value for the 
atmospheric extinction for each night. Since April 1999, the 
measurement is compl eted in the r'-band at an effective 
wavelength of 625 nm (|CMTll2007t) . For a dust-free night, 
the extinction for the passband r' was determined to be 

0. 09 mag. The atmospheric transmission can be calcu- 
lated with the following formula, where too is the extinction 
in a dust-free night and m the observed one: 

L = 1 o r ^£? (i) 

While the Cherenkov spectrum extends between 280 nm 
and 700 nm and peaks at about 320 nm, the optical mea- 
surements are above 600 nm. Since the effect of SAL is dom- 
inated by Mie scattering, which is independent of wave- 
length, the measured absorption can be used to complete 
the correction. 

Since the values of the CMT are provided at zenith, 

1. e. airmass X = 1, they must be corrected for the higher 
airmass when using them for data of higher zenith distance. 



For the IACT data presented here, the zenith distance is 
between 18° and 35° (1.05 < X < 1.22), which corresponds 
to a difference for the atmospheric transmission of less than 
5 % for the range of extinction values shown here. 

3.2. KVA - Optical Data 

The KVA telescop43 consists of two telescopes with 60 cm 
and 35 cm apertures. The latter is used for photometric 
observations of TeV blazars in the R-band (effective wave- 
length 640 nm). The flux of the blazar is determined by 
measuring the counts inside a circular aperture and com- 
paring this to several calibrated comparison stars in the 
field of the blazar, measured with the same aperture ra- 
dius. In this way, any influence of atmospheric extinction is 
effectively eliminated. However, the comparison stars can 
also be used to estimate the atmospheric extinction dur- 
ing the night. In principle, the extinction could be deter- 
mined by fitting a straight line to an airmass-magnitude 
plot; the photometric accuracy of the comparison star data 
in the literature (0.02mag-0.05mag) is, however, insuffi- 
cient for this. Only stars observed at low airmass (X<1.2) 
were instead selected and, for each star, the zero point 
value zq — Rinstr — Real was determined individually, where 
Rmstr was the measured magnitude (scaled to an exposure 
time of one second) and R ca / was the calibrated R-band 
magnitude. Averaging z$ for each night, a value of 5.16 was 
found to provide a good approximation of a photometric 
night. During nights affected by the SAL, zo was higher. 
The transmission relative to a dust-free night was com- 
puted by Eq. [1] where too was 5.16 and to was taken to 
be the average zq for the night. 

3.3. MAGIC Data 

In addition to optical measurements, the atmospheric 
transmission can be determined directly from the IACT 
data. 

A simple way is to use the absolute light yield measure 
obtain ed from muon images as described in iGoebel et alj 
(2005). This method is based on the fact that the spectral 
shape and flux of muons at the telescope are well known. 
For a well-determined muon sample, only showers where the 
telescope is inside the Cherenkov light cone are selected. As 
these showers produce ring-like images in the camera, the 
background can be completely suppressed. Comparing the 
measured light yield distribution to that of a corresponding 
simulated sample, a ratio is determined that can be used 
for an absolute light calibration to correct for changes in 
the light efficiency of the telescope, e.g. mirror degrada- 
tion. Since this correction is completed on a monthly basis, 
short time variations in the atmosphere are not taken into 
account. Consequently, the effect of SAL is still evident in 
the measured muons. Having an impact parameter smaller 
than the diameter of the mirror (17 m), the selected muons 
produce their light between 700 m above the telescope and 
the ground, i.e. between 2900 m a.s.l. and 2200 m a.s.l. 
in case of the MAGIC telescope. To account for the fact 
that this light passes only part of the SAL, the determined 
ratio is multiplied by a factor of ten to estimate the ab- 
sorption, assuming an effective production height of 350 m 
for Cherenkov light. Without any additional measurements, 

http:/ /tur3. tur.iac.es/ 
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1. In the upper plot, the num- 
produc ed photons (genera ted 
with CORSIKA (|Knapp fc Hecldl2004ft ') is 
shown versus height (airmass in g/cm 2 ) for 
monoenergetic showers with energies be- 
tween 50 GeV and 500 GeV in different 
colors. Below, the integral fraction of pro- 
duced photons is plotted versus height for 
the same showers. In addition, the inte- 
gral fraction of light at to 5.5 km, indicat- 
ing the upper border of the Saharan Air 
Layer, is listed for each graph. 



the layering of the SAL is unknown. Accordingly, the results 
from the muons are not very precise and should therefore 
not be used for correction, but can serve as an independent 
crosscheck and first-order approximation. 

The absorption can instead be determined directly from 
the images of the showers themselves by comparing the 
SIZE distributions of the affected data with those of a dust- 
free night. The transmission is calculated from the shift in 
the peak of these distributions, where p is the position of 
the peak either for a dust-free night (j> dust free) or for a 
night affected by the SAL (p a ffected)- 

I Pdustfree Paffected 

-?0 Pdustfr 

For each night of affected data, one value 4- was deter- 



mined. Since the observations used have almost the same 
zenith angle distributions and the zenith distance is only 
between 18° and 35°, the dependence on the zenith angle 
can be neglected. 

A similar method was studied for the Whipple telescope 
(|Lebohec fc Holder! [20031 ) using the measured cosmic ray 
background for a relative calibration of the detector re- 
sponse. 

4. Correction Method 

In the analysis of MAGIC, there are several levels at which 
possible corrections can be applied: 

1. the factor in the absolute light calibration can be 
adapted to correct for the additional absorption. 

2. the SIZE of the shower can be increased according to 
the absorption 



3. the correction can be included in the energy reconstruc- 
tion. 

The last two possibilities obviously do not correct for 
changes in the shower morphology. Losses of light at the 
shower margin might change the shape of the image after 
cleaning. Consequently, common background-suppression 
cuts, normally applied to all data, will give less accurate 
results, since the morphology of the showers is different 
from the one assumed by the cuts. Inter-night changes in 
the atmospheric absorption complicate adapting the cuts 
to the changing shower shape. In contrast, a correction ap- 
plied before cleaning, i.e. to the raw signal of all pixels, also 
correctly recovers the image shape and parameters accord- 
ingly. 

For this approach, an additional factor is applied in 
the conversion from FADC counts to photo-electrons (phe). 
The correction is applied according to the following for- 
mula: 

signal[phc] = signal[FADCcounts] • c ca i ■ y- (3) 

Applying this correction, obviously non-detected show- 
ers cannot be recovered, but the SIZE distributions are not 
affected significantly for an atmospheric absorption up to 
40%. For simulated data (including absorption), the SIZE 
distributions agree well after applying the correction for an 
absorption smaller than 40 %. In cases of higher absorption, 
a more detailed Monte Carlo study is required, and a more 
complex correction might be considered. 

Furthermore, the fluctuations in signals, such as those 
from night-sky background light, increase. Comparing cor- 
rected dark-time observations with data taken in strong 
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moonlight (|Albert et al.l l2007h . the fluctuations after cor- 
rection are found to remain well below those for extreme 
light conditions. By including arrival-time i nformation in 
the i mage cleaning, it could also be shown (Tcs caro et al.l 
12007ft that even data with strong fluctuations are properly 
handled. Hence, no further correction is necessary to ac- 
count for this effect if the correction is applied to data 
taken under normal light conditions. Consequently, the cal- 
culated image parameters are not skewed by the increase 
in the amplitude of the fluctuations and therefore do not 
influence the efficiencies obtained from simulated showers. 



5. Results 

To check the method, a source of constant flux is ide- 
ally used. However, no affected MAGIC data of a steady 
source are available. Data affected by the SAL are diffi- 
cult to schedule, since these rare weather conditions are 
hardly predictable and difficult to identify for the support 
astronomers. 

In July 2006, the MAGIC telescope observed the Active 
Galactic Nucleus PG 1553+113 contemporaneously with 
the X-ray satellite Suzaku, the IACT array H.E.S.S. and 
the optical telescope KVA during a MWL campaign. Due 
to the SAL, the absorption for the single nights from July 
15 th to July 28 th varied between 5 % and 30 % according to 
the SIZE measurement. Therefore, the discussed correction 
was applied to the MAGIC data. 



values from the CMT. For the shown distributions, special 
cuts were applied to remove non-physical and muon events. 
Each curve corresponds to data of one night. 

Before the correction, the SIZE distributions spread be- 
tween 2.8 a. u. and 7a.u. at 200 phc. After the correction, 
the spread was reduced to the interval between 4.7 a. u. and 
7a.u. apart from one outlier, corresponding to one night of 
variable extinction (red dotted line, 22.7.2006). Apart from 
this and a second curve (black dotted line, 17.7.2006) with 
a deviation at around 300 phe, the SIZE distributions agree 
well after the correction for SIZE > 200 phc, i.e. in the range 
used in the analysis. 

5.3. Signal Reconstruction 

The cuts for the background suppression mentioned in 
Sect. I2.l1 include an intrinsic lower cut in SIZE, i.e. events 
with small SIZE are removed. Consequently, events with re- 
duced SIZEs due to atmospheric absorption are cut away, 
resulting in a lower significance of the signal if the same 
standard cuts are applied. 

For PG 1553+113, a significant number of events are ex- 
pected at low energies due to the measured steep spectrum 
of the source. The expected effect of reduced significance 
due to the SAL is observed for data taken in July 2006: 
Without correction, a signal is obtained with 3.2 er signif- 
icance, improving to 5.0 er with correction, using identical 
cuts for the background suppression in both cases. 



5.1. Measurements of the Atmospheric Transmission 

With the atmospheric transmission determined from the 
presented measurements, comparable results were achieved 
for the eight nights of observation of PG 1553+113 in 
July 2006 (Fig. [^. From the MAGIC data, the transmis- 
sion was determined only for the nights in which data of 
PG 1553+113 were available. During some nights, the at- 
mospheric transmission was changing, as seen both by the 
CMT and KVA. The KVA reports considerable changes 
even within a few minutes. Consequently, the different mea- 
surements agree less than for nights with stable extinction, 
since they were probably acquired not precisely simulta- 
neously. Unstable conditions during the observations wors- 
ened the accuracy of the correction. The remaining differ- 
ences between the values from the KVA and the CMT can 
be explained by the fact that the intrinsic errors in magni- 
tude for the comparison stars are larger in the case of KVA. 
Deviations of the SIZE may come from the limited accuracy 
of the method. For the muon measurements, the systematic 
uncertainties in the values are larger due to the strong de- 
pendence on the vertical distribution of the sand-dust in 
the SAL. Apart from small deviations and the differences 
in the nights with variable extinction, the measurements 
agree well. 

Since the CMT measurements are the most accurate, 
they were used for the calculation of the correction factors. 

5.2. SIZE Distributions 

As explained in Sect. [2~T1 the atmospheric extinction influ- 
ences mainly the SIZE of the showers. Figure [3] shows the 
SIZE distributions for the range used for the analysis be- 
fore (left) and after (right) applying the correction with the 



5.4. Flux Determination 

Applying the correction as explained in Sect. 21 the 
8.5 hours of MAGIC data of PG 1553+113 from July 2006 
were analyzed. Fitting a power law to the differential spec- 
trum, aflux of (1.4±0.3)-10~ 6 phTeV- 1 s- 1 m- 2 at 200 GeV 
and a spectral index of —4.1 ± 0.3 were determined. A 
daily light curve shows a flux consistent with a constant 
fl ux within the error s. The results were shown in detail 
in I Albert et all (|2008[ ) and were used for MWL studies by 
iReimer et all (|2008h . 

Correcting the flux for the effect of the SAL, the sys- 
tematic uncertainty in the method of 5 % results in an ad- 
ditional error of 0.2 • 10~ 6 ph TeV~ 1 s~ 1 m~ 2 in the absolute 
flux. 



6. Conclusion and Outlook 

The different measurements of the atmospheric transmis- 
sion from stellar photometry (KVA, CMT) and IACT data 
(MAGIC) agree well for nights with stable extinction. Since 
more than 90% of the shower light originates above the 
SAL, the measurements can be used directly to correct the 
IACT data. By applying a correction factor in the cali- 
bration, the data of PG 1553+113 taken during the MWL 
campaign in July 2006 for example were corrected success- 
fully. This shows that the loss of observation time due to 
atmospheric absorption can be reduced. 

For additional checks, affected data of a steady source 
are required. When such data become available within 
scheduled MAGIC observations, the presented method can 
be tested further. 

The correction should not be applied blindly to data 
affected by atmospheric absorption of another kind. For 
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Fig. 2. Atmospheric Transmission for the nights between July 16th and 28th 2006 from the measurements of the CMT 
(red dots), the KVA (blue full squares) and MAGIC. From the latter, not only the SIZE measurement (black open 
squares) is shown, but also the values derived from muon rings (green circles) where the value is scaled by a factor ten to 
account for the smaller absorption for muons. Nights with unstable extinction are marked with gray shaded background. 




Fig. 3. SIZE distributions between 200 phe and 2000 phe (i.e. approximately from 100 GeV to ITeV) before (left) and 
after (right) correcting for the effect of SAL. Each line corresponds to one night and is scaled by the observation time. 
The color and style of each line is different to distinguish the nights from each other. 



those, the vertical distribution of the atmospheric condi- 
tions should be studied e.g. with a Lidar (light detection 
and ranging), i.e. a device shooting a laser beam into the at- 
mosphere and measuring the backscattered light. The pos- 
sible use of information provided by a Lidar, however, still 
has to be investigated. 

In cases of strong absorption, i.e. above 40%, the 
method is limited by the fluctuations in the night-sky back- 
ground light and the effect cannot be corrected completely. 
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